Black hole interacting with matter as a simple dynamical system 



Petr Hajicek 
Institute for Theoretical Physics 

University of Bern 
Sidlerstrasse 5, CH-3012 Bern 
Switzerland 
(February 7, 2008) 



Abstract 



Recently, a variational principle has been derived from Einstein-Hilbert and a 
matter Lagrangian for the spherically symmetric system of a dust shell and a 
black hole. The so-called physical region of the phase space, which contains all 
physically meaningful states of the system defined by the variational principle, 
is specified; it has a complicated boundary. The principle is then transformed 
to new variables that remove some problems of the original formalism: the 
whole phase space is covered (in particular, the variables are regular at all 
horizons), the constraint has a polynomial form, and the constraint equation 
is uniquely solvable for two of the three conserved momenta. The solutions 
for the momenta are written down explicitly. The symmetry group of the 
system is studied. The equations of motion are derived from the transformed 
principle and are shown to be equivalent to the previous ones. Some lower- 
dimensional systems are constructed by exclusion of cyclic variables, and some 
of their properties are found. 
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I. INTRODUCTION 



The spherically symmetric gravitating shell has been used as a simplified quantum model 
at many occasions. One can mention the quantum effects of domain walls in the early 
Universe [p], quantum aspects of gravitational collapse @, quantum theory of black holes 
or Hawking evaporation of black holes ||]. 

The equations of motion for the shells has been derived by Dautcourt || and transformed 
to a geometric (i.e. "gauge invariant") form by Israel |J. However, for a quantum theory, 
one needs a variational principle rather than dynamical equations. Such variational principle 
has as yet been just guessed from the dynamical equations, Refs. |7|] and 0, or from some 
intermediate variational principle, Refs. [|j] and §J . Quite a number of different Hamiltonians 
have resulted and the corresponding quantum theories have not been unitarily equivalent. 
There are two sources of the ambiguity. 

The first source of ambiguity is the invariance of the general relativity with respect to 
any coordinate transformations on one hand, and the property of Hamiltonians to generate 
dynamics with respect to a particular time coordinate on the other. For example, the choice 
of the proper time along the shell as such a coordinate leads to the coshp Hamiltonian, 
Refs. and [ lOfl , the Schwarzschild time coordinate inside the shell to the square-root 
Hamiltonian, Ref. 0, and the Schwarzschild time coordinate outside the shell to a merely 
implicitly determined Hamiltonian of Refs. [[| and [[|. An attempt to work out a gauge 
invariant theory for at least some class of the guessed Hamiltonians is Ref. [^]; based on 
a technical assumption that the super-Hamiltonian is quadratic in the momenta, it gives a 
unique action principle. 

The second source of ambiguity is the fact that equations of motion do not determine 
the corresponding variational principle uniquely in general. This ambiguity can be removed 
by a direct derivation of the variational principle for the shell from the Einstein-Hilbert and 

they all lead 



|g] and 
and 
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13|| , the reduction 



shell matter Lagrangian. Three such derivations exist: Refs. [|Tl] 
to the same variational principle for the same system. In Refs. 
of the second order formalism to the spherically symmetric case is done first, followed by a 
reduction to dynamical variables of the shell alone and by the transformation to the second 
order formalism. In Ref. |ll[j, a first and second order formalism for a general shell and 
gravitational field is derived (no symmetry); this second order formalism is reduced by the 
spherical symmetry and to shell variables in Ref. 
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We are then left with the following dilemma. The super-Hamiltonian of Ref. || is simple 
and amenable to the existing gauge invariant quantization methods (like e.g., Refs. []T7| and 
The super-Hamiltonian of Refs 



121, 14 and IT 



follows from the Einstein-Hilbert 
and shell-matter Lagrangian, but it is extremely complicated: it contains nested square 
roots and hyperbolic functions, and for a good measure, it is formulated in coordinates that 
diverge on horizons. It would be difficult to quantize. 

A natural question then arise: what is the relation between the guessed and simple 
variational principle of Ref. || and the derived but complicated one of Refs. [12], \TA\ and 
|13|| ? If the simple principle were equivalent to the complicated one, we could use the 
simple principle and forget the complicated one. However, the answer turns out to be rather 
surprising (see Ref. |TtJ): the two descriptions are only locally, but not globally equivalent. 
The local equivalence explains, why the equations of motion are the same. For the study 
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of the problem, a geometric approach of Ref. has been applied. All gauge invariant 
properties of a constrained system are encoded in two objects: the constraint manifold T 
and the presymplectic form Q on T. The equivalence of two such systems, (r^f^) and 
^2,^2) is then a well-defined mathematical concept independent of a choice of extended 
phase space, a choice of constraint functions or a choice of gauge. 

In this situation, it seems natural to look for a transformation of the variational principle 
of Refs. [12], [13 1 and [jnj to a better set of coordinates, and this will be the main topic 



of the present paper. The method will again be based on the gauge invariant description 
(T, Q) of reparametrization invariant systems. We shall perform the transformation in three 
steps; to motivate the steps, we take into account the symmetry of the system, the structure 
of the so-called Cartan form B, which is defined by d® := Q, and the topology of the 
space (r, f2). The topological problem involved here is the existence of transversal surfaces 
in T; such surfaces are nowhere tangential to the dynamical trajectories. Functions whose 
levels are transversal are not only helpful, if one looks for nice coordinates, because they 
simplify B. Quite generally, they have to do with the existence of Hamiltonian and with the 
possibility to give the quantum dynamics the form of Schrodinger equation (cf. Ref. [PHI). 
We shall, therefore, also look for transversal surfaces in a more systematic way. Finally, 
we shall find that the structure of the phase space of the spherically symmetric shell is not 
simple. Bnly a proper subset, which we call physical region, contains physical states of the 
shell. Its nature can be roughly specified as follows. All points inside it correspond to the 
system consisting of a shell of positive mass and radius, interacting with a black hole. The 
boundary of the region is complicated; it contains shells with zero radius, shells with zero 
rest mass (which cannot, however, be regarded as shells made of light-like matter), as well 
as self-gravitating isolated shells with positive radius and mass — these are the points at the 
boundary that are physically meaningful. The points outside the region that describe shells 
with positive radius and mass correspond to systems consisting of the shell and a negative 
mass source. A conclusion seems to be that the system will be difficult to quantize even if 
the new variables has made it algebraically simple. Bne problem will be to keep the spectra 
of observables contained within the physical region. Another problem will be to construct 
a unitary dynamics, because the classical dynamics breaks down at the boundary of the 
physical region. 

The plan of the paper is as follows. In Sec. [TT], we collect the relevant results of the 
previous papers, mainly Refs. [14] and [|H| . The variational principle that we describe allows 
also for an internal degree of freedom of the shell; it depends, therefore, of two additional 
variables in comparison with Refs. [|1^] and []14" |: the proper time along the shell trajectory 



and the rest mass of the shell. We determine the constraint manifold T and the Cartan form 
B that follow from the variational principle. We study the physical region. 

In Sec. H|, we introduce three new functions on the phase space that are regular at 
the horizons and that will replace the momentum P of the shell and the two Schwarzschild 
times T ± that diverge there. We describe their physical and geometrical meaning, using 
also some results of previous papers. In Sec. [TV], we specify a new set of coordinates on 
T, and show that the transformation from the old to the new coordinates is differentiable 
and invertible. We also find the ranges of these coordinates in the physical region. Sec. [V] 
collects and extends the results about symmetries of the system that were obtained in Ref. 
[IJ| and expresses the symmetry transformations in the new coordinates; the form of the 
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transformation simplifies. The new times are cyclic coordinates as the old have been. In Sec. 
VI , we transform the Cartan form into new coordinates. The square roots and hyperbolic 
functions disappear, and an equivalent form becomes regular at the horizons so that it can 
be smoothly matched accross them. However, the Cartan form is still complicated and this 
motivates further transformations. The final result thereof can be described as follows. The 
extended phase space is (R 8 ,f2). The natural coordinates of R 8 form a Darboux chart for 
Q. The equation of the constraint surface T in R 8 is polynomial in all momenta and the 
presymplectic form dQ is the pull-back of Cl to T. 

In Sec. |V11] , we study the polynomial constraint. We show that it defines two disjoint 
submanifolds in R 8 . One, denoted by r 6 , is six- dimensional and lies well outside the physical 
region. The other, denoted by Tj, is seven-dimensional and intersects the physical region. In 
the intersection of T-? with the physical region, we find two different foliations by transversal 
surfaces. They are not globally transversal, however, because none of them is intersected by 
all dynamical trajectories. The problem is that each dynamical trajectory starts or finishes 
at the (R = 0)-subset of the boundary of the physical region. The two different Hamiltonians 
coresponding to the two foliation are written down explicitly; they contain second and third 
roots. 

In Sec. |V11J| , we repeat the study for some important lower- dimensional cases; they can 
be obtained by exclusion of cyclic coordinates from the original system. In particular, a 
(dynamical) black hole and the shell without its internal degree of freedom, the shell in the 
field of a fixed black hole external field, or the isolated self-gravitating shell (with the flat 
spacetime inside) are considered. 



II. DESCRIPTION OF THE SYSTEM 

In this section, we shall collect some results of Refs. ||, |TJj] and [[HJ so that the paper 
becomes relatively self-consistent. We shall also describe and study the physical region. 

A spherically symmetric thin-shell spacetime solution of Einstein equations can be con- 
structed as follows. Consider two Schwarzschild spacetimes M.\ and M.2 with Schwarzschild 
masses E\ and Ei. Let Si be a timelike hypersurface in M.\ and £ 2 be one in .M 2 . Let 
Ei divide M.\ into two subspacetimes, M.\+ and -Mi_, and similarly E 2 divide M.2 into 
M.2+ and M.2-- As everything is spherically symmetric, all spacetimes are effectively two- 
dimensional; we chose fixed time and space orientation in the two-dimensional Schwarzschild 
spacetimes so that future and past as well as right and left are unambiguous; let then .M 2 + 
and be right with respect to M.2- and Mi-. Let Ei and E 2 be isometric; then the 

spacetime Mi- can be pasted together with the spacetime M2+ along the boundaries E x 
and E 2 . The result is a shell spacetime M s . Given a shell spacetime, we shall leave out the 
indices 1 and 2, having right (left) energy E + {E_), shell trajectory E, and the right (left) 
subspacetime M+ (M-). Thus, M+ = M 2 C M s and M- =Mi(~) M s . In Refs. (TTJ, [T5| 
and ||, the Schwarzschild masses E e , e — ±1 were assumed to be positive. 

Each subspacetime M e has the metric 

ds 2 = -F e (R)(dT e f + F e - 1 {R)dR 2 , (1) 

where 



4 



F e {R) : = 1 - 




and it is split up by its horizons into four quadrants: we denote by Qi that which is adjacent 
to the right infinity, Qu to the left infinity, Qui to the future singularity and Qiy to the 
past singularity. Suitable notation (introduced in Ref. enables us to write all formulas 
in a form valid in any quadrant: we define four sign functions a e and b t distinguishing the 
quandrants of A4 e : a t := sgnF e , and b t equals to +1 (—1) in the past (future) of the event 
horizon in the spacetime (A4 e ,g t ). We also include the dust internal degree of freedom, M, 
and the proper time T along the shell trajectory using the results of Ref. [|l^]. Then the 
total action for the dust and gravity can be rewritten in the form 

S E [T, M; T ± , E±\ R, P;u] = J dt(PR - E + f + + P_T~ + MT- uC), (2) 

where v is a Lagrange multiplier, 



C = M - P^F + + F_ - 2a_M_\/|F + F_|sh a+a _^, (3) 

is the super-Hamiltonian, sh a x := (e x + ae~ x )/2 for any a = ±1 and x G (—00,00). The 
momentum P conjugate to the radial coordinate R can be defined as follows (for more detail, 
see Ref. ||14]| ). Let (n e ,m e ) be the orthonormal dyad at a shell point in the subspacetime 
(Ai e ,g t ), the vector n e being tangential to the shell and future oriented, m e being right 
oriented; we call it shell dyad. Let further the so-called Schwarzs child dyad (n Se ,ms e ) be 
defined at each point of the shell in the spacetime (A4 e ,g e ) by the requirement that n St be 
timelike future oriented, ms € be right oriented and that one of the vectors be tangential 
to the (R = const)-curves. Clearly, in each quadrant, we have a different formula for the 
components of the Schwarzschild dyad; with respect to the Schwarzschild coordinates, 

for F e > 0, n Se = (b e /^F e , 0), m Se = (0, h^E); 

forF e <0, n S€ = {0,b £ yJ\K\), rn Se = {-bJy/\F\\, 0). 

Then, P e /R is the hyperbolic angle between the Schwarzschild and shell dyads: 

n e = rise coshP e / R + ms t sinhP e / R, 
m e = rise sinhP e / R + m St coshP e / R. 

Finally, P := [P], where we use the common short-hand for the jump [A] := A + — A_ of a 
quantity A across the shell. 

The extended phase space of the system is eight-dimensional, split up into 16 disjoint 
sectors, each sector being a pair of quadrants, one chosen from the left and one from the 
right subspacetimes. Each of these sectors can be covered with the coordinates P, R, E + , 
T + , T~ , M and T. The constraint surface that we call T is defined by the constraint 
equation C = in the extended phase space. It does not intersects all 16 sectors, but it is 



5 



split up into at most (if E- > M) non-empty intersections of the sectors with it.[] These can 
be covered by the system of seven coordinates P, R, E + , T + , E_, T~ and T. The following 
coordinates are assumed to have non-trivial ranges: 

R > 0, E + > 0, E_ > 0, M > 0. (4) 

The pull back of the Liouville form of the action (Q) to the constraint surface T in these 
coordinates is simply 

= PdR - E + dT + + E^dT- + MdT, (5) 
where M is now a function of the seven coordinates defined by 



M = R ] Jf + + E_- 2a_b + b^\F + F_\sh a+a _^. (6) 

The form (J|) is called Cartan form of the system; it contains all information about the 
dynamics and about the Poisson brackets (see, e.g. Refs. |I| or JT5J). 



The present paper starts from these formulas and quantities. The first new observation 
that we make is, that all these formulae remain valid also for the values E + = or EL = 0, 
if the notation and conventions are adapted a little. These cases are interesting because 
they describe an isolated self-gravitating shell (with the flat spacetime inside). Indeed, the 
metric of the flat spacetime is covered by formula ([I]). However, the flat spacetime does not 
contain any horizon and is not split into quadrants. Still, a piece of the flat spacetime can be 
used in four distinct ways in the contruction of the shell spacetime: either it lies to the left 
(E_ = 0) or to the right (E + = 0) of the shell, and, in each case, the R coordinate can either 
increase to the right or to the left (R is always spacelike in the flat spacetime). We can, 
therefore, formally define the E-spacetime for E > as before to be the Kruskal spacetime of 
Schwarzschild mass E with four quadrants, and for E = as consisting of two topologically 
separated quadrants Qi and Qu, each isometric to the flat spacetime, but with oposite time 
and space orientations: in Qj, the radial ("Schwarzschild") coordinate R increases to the 
right and the "Schwarzschild" time T to the future; in Qu, R increases to the left and T 
to the past. It is amusing to observe that the validity of the formulae @ and (§) can then 
be extended to vanishing E's, if we just use the values a e = 1 everywhere, b e = +1 in Qj 
and b t = —1 in Qu. This works because the components of the Schwarzschild dyad that 
are crucial for the derivation of Eq. (|3p in Ref. [fTl] retain their form in the quadrants Qi 



and Qu of the .E-spacetime for E = 0, and there are no quadrants Qui and Qiy. Similarly, 
the sign relation (14) of Ref. ||, which is basic for the form of the dynamical equations, 
remains valid with our orientation and quadrant conventions. Observe that the presence of 
two topologically separated subspacetimes is not observable and does not imply any physical 



1 Observe that we admit unphysical sectors, i.e. those, in which the constraint equation C = has 
no solutions. This simplifies the boundary of the extended phase space a great deal and does not 
lead to any problem: it is the very purpose of extended phase spaces to enclose also non-physical 
poits, if this leads to any simplification. 
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assertion; similarly, the orientation of coordinates is a mere coordinate convention without 
any physical meaning. 

A very important but rather embarrasing observation is now that a formally impeccable 
shell spacetime can be pasted together from Schwarzschild spacetimes of any mass, even a 
negative one. The metric is still given by Eq. (|1|) and the spacetime has a similar global 
structure as the flat one: there are no horizons and there are two possible orientations of how 
it can be built in. Again, all equations remain valid, if we accept the same sign, quadrant and 
orientation conventions as for the flat spacetime and define the (two-quadrant) i?-spacetime 
for E < analogously to what we have done for E = 0. The corresponding shell dynamics 
is then again given by the action (fj). 

Let us quickly discuss the dynamical equations following from the action (§) under these 
conventions. There are three conservation laws, 

E + = E_ = M = 0, (7) 

resulting from varying the action with respect to the cyclic coordinates T + , T~ and T. The 
variations with respect of the conserved momenta and P with the subsequent simplification 
by the constraint (H) yield: 

j = u, (8) 

(9) 




= ua _ b+b __^\ F+F _\ sh _ a+a _- (10) 

The equation determining P can be obtained either by varying the action with respect to R, 
or by differentiating the constraint equation C = with respect to t followed by substitution 
for all other t-derivatives from Eqs. (|7|)-(]T0|). Let us rewrite the constraint equation in the 
form 

MR r p 

+ -r-.(F+ + F\ = a-b+b-— J\F+F J sh a + a — . (11) 

2i? 2M v+ ; + M VI + 1 + R v ; 

Eqs. (pD-JTID represent a complete system of dynamical equations for the shell. 

We derive some important consequences from the dynamical equations. Eqs. fl8D— ([LTD 
imply that 

„ (dT e \ 2 1 (dR\ 2 

confirming that T is the proper time along shell trajectories. Eqs. @, (|J) and ( |TTD deliver 

the time equation: 

_ alT e M E + — E_ 

~dT = " e 25 + -^TS — ■ < 12 > 

From Eqs. flTU| ) and (0), the radial equation follows, 
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where 



V(R) 



M 2 

4R2 



E + + E_ (E+ - E_ 



R 



M 5 



+ 1 



(14) 



and uj = ±1 is a suitable sign (see Eq. (2) of Ref. fll5"f). The square of the radial equation 
can be rewritten in one of the two forms 



'dR s 
df. 



according as e = ±1, where 



J, 



M 



J 



E,-E_ 



e ' 2R 6 ~ 



M 



(cf. Ref. ||). It follows immediately from this definition of J e that 

J+ + J- 

and this can be transformed with the help of Eq. ( |i~5|) to 



M 



sgnJ 4 



\ 



r. ( dR Y 



M 



'dRY _ _ 
dJ J ' R 



However, comparing Eqs. (|i~2"D and (0), we observe that 

sgn J e = — e sgn F e 



dJ 



The sign of the expression F e (dT e /d~T) deserves a name; let us call it r e (cf. 
obtain finally: 



'dR s 
df. 



+ r_ 



\ 



dR^ 
dJ , 



M 



(15) 



(16) 



Then, we 



(17) 



This is the so-called Israel equation for dust shell (cf. Eqs. (14), (21) and (22) of Ref. ||). 

Eqs. ([L2f)-(ll7D are valid for all values of Schwarzschild masses E + £ (— oo, +oo) and 
E_ £ (— oo, +oo). Some features of the unphysical (negative E £ ) cases are described by the 
following theorem. 

Theorem 1 Let E e be non-positive and not larger than E^ e for some value of e. Then the 
timelike (connected) center (R = 0) curve lies in the subspacetime Ai e on the e-side of the 
shell. 
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To prove the theorem, we observe that E e < E_ € implies 




Thus, the left-hand side of Eq. (|TT|) can only be positive, if r e = — e. The assumption 
E € < of the theorem implies that A4 e is the two-quadrant E'-spacetime with r e > in the 
quadrant Qi and r e < in the quadrant Qn; this follows from the definition of r e and from 
the positivity of F e (in fact, for E e < 0, we have F t > 1). Hence, to the right of the 

shell, there is a part of Qu with the center to the right, and to the left of the shell (A4-), 
there must be Qi with the center to the left, Q. E. D. 

An interpretation of the theorem is that the shell spacetime containing one or two 
negative-mass subspacetimes is unphysical: there must be at least one negative mass source 
somewhere outside the shell. The theorem also says that a shell spacetime constructed 
from one flat and one positive-mass subspacetimes must have a regular center inside the flat 
subspacetime. 

To summarize: The action (^) generates a regular dynamics in the enlarged phase space 
defined just by 

M>0, R>0. (18) 

For E + > and E_ > 0, it describes a shell interacting with a black hole of mass E_ (or 
E + , depending on from where we observe the shell). For E + > and E- = (or E + = 
and E- > 0), it describes a self-gravitating shell with flat spacetime inside. If any of E + 
and E_ are negative, it describes the shell interacting with a negative mass source. Only 
the points with 

E + > 0, E_ > (19) 

are physically sensible; the subset specified by the inequalities (jig) and (|i"9f) is the physical 
region of the phase space. There seems to be nothing about the variational principle (0) 
that would help us to enforce the validity of the Eqs. (|i~9|) in any natural, automatic way. 



III. RADIAL AND KRUSKAL MOMENTA 

From the definition of the functions and P as given in the previous section, it follows 
that they are singular at the horizons: T e at R = 2E e and P at R — 2E_ as well as 
at R = 2E + . In this section, we introduce three functions that are well-defined on the 
whole constraint surface to replace T ± and P. Another replacement of and P with 



similar properties has been tried in Ref. [14|. There, the functions were constructed from 



the Kruskal coordinates; however, the Kruskal coordinates make no sense for flat spacetime 
and so the important case E- = could not be incorporated. In this subsection, we remove 
this problem. 

Let us define the new functions q, [Ti] and T\ by 
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f x :-- 



T+ + T- 



[Ti] :=T + -T- + 2E+\n 

and the auxiliary functions u t and v e by 

M 



sn_ a 


p 


1 I 


u + 


+ E_ In 


v+ 




u + 


- 2E_ 







u. 



(20) 
(21) 

(22) 



q__[E] 
M M 



2R J 



v f : = - — 



I + IS 
M M 



M 

'2R' 



(23) 



in Eq. Q2"3]), the functions g and M are given by Eqs. (^j) and := — 

The meaning of the momentum q can be seen, if it is expressed by means of velocities 
along dynamical trajectories. Comparing Eqs. (W), ([II]) and (PD|), we find that 
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M 



dR 



(24) 



hence, g is a kind of radial momentum, and it is regular everywhere along dynamical tra- 
jectories. 

The meaning of the functions T\ and [Ti] can be inferred from that of Tf, which are 
defined by 



T 1 -=f 1 + -[T 1 ]. 



Eqs. 



and 



lead then to 



TI = T e + 2E t In 



for each e = ±1. Let us limit ourselves to the physical region and distinguish two cases. 
E € = . Then, T{ = T € in both quadrants of the -E-spacetime, so b € T[ is the time coordinate 
of the inertial system of the Minkowski spacetime Ai e in which the center of mass of the 
shell is in rest. E e > . Then, it is possible to introduce Kruskal coordinates U e and V e in 
the Kruskal spacetime A4 e by 



R = 2E € k{-U £ V € ), 



T £ = 2E f In 



Ur 



(25) 



and by the requirement that both functions U e and V t increse to the future, where k is the 
well-known monotonous function defined on the interval (—1, oo) by its inverse, 

/t _1 (x) = (x - l)e x . 

At each point of the shell trajectory U e = U e (t) and V 6 = V e (t), the so-called Kruskal 
momentum, P^, was defined in Ref. by 
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P* -^ln^ 
Fk 2 ln dU e 



In Ref. ||15|| , four further functions u e and v e were defined by 



(26) 



where K f 



-U e V\ 



The following relation between the pairs u e , v e and R, E e was shown to hold in Ref. [15 
(Eqs. (32) and (33)): 

[ E \ M r-T? 
i E ] M r-T? 



A comparison with Eqs. (|13|) and (|14D shows that our functions w e and v e defined by Eqs. (p3|) 
coincide with the functions M e and f e of Ref. [|T5 1 , if the dynamical equations are satisfied. 
Eq. ( p6| ) then shows that u € and f e vanish at horizons and that their logarithms diverge 
there. It then also follows for the functions T{ along dynamical trajectories that 



R ' 



(27) 



We can see that the functions are directly determined by the geometry of the dynamical 
trajectory in the spacetime Ai e and are regular everywhere along the trajectory. 



IV. A REGULAR COORDINATE SYSTEM 

In this section, we describe a coordinate system that is regular at the horizons, and at 
the same time remains meaningful for the special values E e = of Schwarzschild masses. 

As such coordinates on T, we suggest the seven functions q, R, M, T, [E], 7\, and pi]. 
Let us show that the Jacobian of the transformation from these coordinates to P, R, E + , 
T + , E_, T~ and T is non-zero. To calculate the determinant, we observe that the four 
columns R, T + , T~ and T of the Jacobian 

d(q,...,[Ti]) 
d(P,...,T) 

contain each at most two non-zero elements; if we expand the determinant along these 
columns, we obtain 

[Ti]) = _ dq_ (dM_ dM_\ dM (dq_ _dq_\ 
d(P, ...J) ' dP \dE+ dE_) dP \dE + dE_)' 1 ' 

For the derivatives involved here, Eqs. (|6|) and ( p0|) yield immediately: 
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dq 



p OM 



— = a_M_^|F+F_|sh a+a _-, Qp = -a-b+b-\/\F + F^\—sh_ a+a _-, 



-R 



dM dM R 



F + + F_ 



dE+ dE_ M 



dq dq 



dE+ dE_ 



i i F + + F - x. 
-a + o + o_ sh. 



F+F_ 

P 



F 



Substituting this in Eq. (|28|) and using again Eq. @, we arrive at the simple result 



d(q,...,[Ti]) 
d(P,...,T) 



M 

R' 



which holds in each sector of T. 

It follows that the transformation is regular within each sector and can be inverted. In 
fact, it is not difficult to find the functions defining the inverse transformation. First, we 
show that 



e = r( q^_[Ei_wr 

' 2 { M 2 M 2 4F 2 



(29) 



in each sector. To this aim, we substitute for q and M from Eqs. (gj) and (|20|) into Eq. fl29| 
and use the identities 



2F X 2 



1 R , 



[EY 
R 2 ' 



(30) 



In this way, we obtain 



2E _ R 2 
~R ~ 1 ~ M 2 



-2(1 - 2F/F) 2 + 2a_M_(l - 2F/F) v /|F + F_|sh Q+a _|J 



Now, the equality follows immediately from Eq. (Q). 

The next non-trivial part of the inverse transformation is the relation 



P=-ln 
2 



U + V- 



(31) 



which again holds in each sector. To show Eq. (pT|), we use Eqs. (|23|) to obtain 



u+v - 



The substitution for q and M gives 



M 2 



9 + r/— = r)(R - 2E) - V a^b + b^\F + F^ p / R , 
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where rj = ±1. After some rearrangement and applying Eqs. fl3~0|), we have 
(q + V^j ~ [E] 2 = -a_ v b + b_R 2 ^\F + F^ p l R (f + + F_ - 2a^b + b + 
Then, Eq. (Of) leads immediately to 



sh a+a _ 



(q + V^j ~ [Ef = -a v b + b^^F~\e^ R , 

and this implies Eq. fl3~l"|). 

The rest of the inverse transformation is easy: the functions R and T are the same in 
both sets of variables, 

E e = E+ € -[E} (32) 

and 

— 6 U 

T e = T 1 + ~[T 1 ]-2E e ln -i ; (33) 

in the last equation, one has to substitute Eq. (|32|) for E e and Eqs. (p3|) for w e and t> e . 

The ranges of the variables q, R, M, T, [i?], 7\, and [Ti] in the physical region are implied 
by the conditions (|T8f) ( |i~9l) and Eq. fl29|) . Let us work them out. Eq. (|29|) together with 
i£ > implies that 

[^] 2 <M 2 + g 2 -^; (34) 

hence, 

M2 + ^>^- (35) 

The two inequalities E + > and E- > are equivalent to [E] 2 < 4E 2 . Eq. ( f29|) implies 
that this inequality is, in turn, equivalent to 

(W - (M 2 + + ^) ) 2 > >/M^) 2 • ( 36 ) 
From Eq. (0), it follows that 



r l2 /\.2 2 M 4 \ M 4 



Eq. fl36|) is, therefore, equivalent to 



[^f< (v / M 2 +^ 2 -^) • (37) 

Eq. ( |35| ) shows that the inequality fl3~T|) is not weaker than (p4|). Hence, all information is 
contained in the following inequality: 



2 



\[E\\<JW + q*~. (38) 

Inequality (|38"D together with M > and R > define the ranges of the variables q, R, M, 
T, Ti, and [Ti] in the physical region. 
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V. SYMMETRY OF THE SYSTEM 



In this section, we collect the results of Ref. [TE\ on the symmetry of the shell action and 
describe the action of the symmetry transformations on the new variables. 

In Ref. [15], we have found that there is a continuous symmetry group generated by an 
arbitrary function of E + and EL. A simple generalization of the argument given in Ref. [TJ[ 



implies that the following finite transformation of the variables E e , u e , v e and P e , which were 
used there, is a symmetry of the system: 

dA 

the other variables being invariant, where A(E + ,E_) is an arbitrary smooth function of 
two variables (the factor 1/4 is introduced for convenience). In addition to this continuous 
infinitely dimensional group, there were two reflections, which we denote here by <j\ and 02. 
o\ is a time reflection defined by 

E € i-> E e , u e + u e h-> —(u e + v e ), u e -u e ^u e -v e , P e ^-P e , 

and o"2 is a left-right reflection defined by 

E + <-> E_, u + <-> V-, m_ <-> v+, P+ «-> — P_. 

These results can be easily expressed in our variables q, R, [E], 7\, E and [Ti] and 
extended to M and T. Consider the function q; Eqs. Q23 ) yield 



g = -M(m £ + u e ), e = ±1. 

It is clear that M is invariant with respect of the whole group; hence, q is invariant with 
respect to the A-transformation for any A(E + , E_), it changes sign by the time reflection 
a 1 and is invariant with respect to Oi- The function R was written in terms of u e and v e in 
Ref. [|15| in the form 

2E e 

R=— — , e = ±l. 

1 + W e f e 

Realising that 4u e v e = (u e + v € ) 2 — (u e — v e ) 2 , we can see that R is an invariant of all above 
transformations. The function [E] transforms non-trivially only by cr 2 , 

[E] i-> —[E], 

and E is an invariant like R. A comparison of our Eq. (E7h with Eq. (43) of Ref. reveals 
that Tl = 4P e , so that its A-transformation is 

Using this formula, one finds easily that 
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f ^~wv |ril " |Ti| -I- (39) 

The time reflection <j\ changes the signs of both times: 

% » -f 1; [7\] ^ -[Ti], 

but the left-right reflection 02 acts as follows: 

f x ^ -T 1; [T x \ » [Tt]. 

From the physical meaning of M and T, it follows that they are invariant with respect to 
the continuous transformation for any A(E + , E_) and of the left-right reflection 02, whereas 
the time reflection a\ must give 

M i-y M, Th-T, 

However, as T is a new cyclic coordinate, the new system has larger continuous symmetry: 
we can extend A to an arbitrary function of the three variables EL, EL, and M, and define 
the action on T as follows: 

T " T + M- < 40 > 



Let us postpone the proof that this extended A-transformation is a symmetry to Sec. VI 



It generates three independent constant shifts of the cyclic variables T\, [TjJ and T that are 
different in different shell spacetimes. 

We can see from these results that our new variables transform particularly simply. This 
was, in fact, the idea that helped to find the variables in the first place. 



VI. TRANSFORMATIONS OF THE CARTAN FORM 



Let us transform the Cartan form G in each sector of T to the variables q, R, M, T, [E], 
Ti and [T x ] and check that it can then be extended smoothly accross the boundaries of the 
sectors. 

If we substitute for P, T + , T~ and M, we obtain 



6 = —Ed[Ti] - \E\dTr + MdT + 2E 2 + d [ In 
u + v. 



u 



+ 



+ In 



U-V + 



d(Ei 



In 



U- 



The logarithms can be rearanged as follows: 



= -Ed[Tx\ - \E\dTx + MdT + d[ ( — + E% ) In 



4 



El )d In 



4 

R 2 



u+ 



2E 2 _d f In 
d(E 2 _). 



^- + El)lu 



U- 



4 
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The total differential of a function (that diverges badly at the horizons) can be left out. We 
obtain 



= -£d[Ti] - [E]dTt + MdJ 



E 2 \d In 



(41) 



where we must substitute for E + , EL, E, u e and v e from Eqs. (p^), (|29D and (|23l). This form 
of G has been obtained in Ref. ]nj (without the matter term); there, however, w e and t> e 
were regarded as independent variables. Observe that Eqs. (0), (p9j ) and (|32"D imply 



2E e 



(42) 



so that w e and v t are not independent functions for E e = 0. 

Now, we can show that the last term in O is regular: we just rewrite it using Eq. 

2E\ 

1 + -^-J (vdu — udv) 





u 






~R 2 i 






\ 








V 






4 V 



(43) 



Thus, our new variables cover all of the constraint hypersurface T as promised. 

Let us express explicitly by means of the variables q, R, M and [E\. Employing Eqs. 
(E3|), we have 



vdu — udv 
[vdu — udv] 



[E] 



2 (±. dq + .L dR -j!L M 



Then, applying the well-known identity [AB] = A[B] + B[A] that holds for jumps and 
averages of any two functions A and B, and after using Eq. (^), we arrive at 




E 2 )d In 



+ R + 



Rtf_ 
2M 2 



u 

V 

3R[E] 
2M 2 



Rq\E\ 
M 2 



d[E] 



2Rq Rq 3 



Rq[Ef 



M M 3 

M 2 \ , / g 3 q[E} 2 

8R dq+ l q+ 2W-m*- 




dM 



dR. 



(44) 



This, together with Eq. (f|l]) gives G as a function of new variables. It is a complicated one, 
so we have to transform the coordinates further in order to get a simple expression. 

The following question will give some direction to this search for simplicity: is the variable 
[Ti] suitable for the role of time? More concretely, the time levels in the constraint surface 
should be transversal (this is shown in Ref. ||10|| ). A surface is transversal, if dynamical 
trajectories intersect it transversally and only once. We can see if the surface [Ti] = const 
has this property as follows. The direction of motion on the constraint surface coincides 
with the direction of degeneration of the presymplectic form dQ; hence, the pull-back dOf 
of dQ to any transversal surface / = const must be non-degenerate. If the surface is 2n- 
dimensional, then the 2n-form dOf A ... A dOf must be everywhere non-zero. In our case, 
n = 3. Let us calculate dO^] A d©[Ti] A dO^]- Eqs. ( f4"l"D and (|4"4j) imply that 
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dQ [Tl] = -d[E\ A dTi + dM A dT 

+ A x dR A dq + A 2 dR A d[E] + A 3 dR A dM + A 4 dg A d[£] + A 5 dq A dM + A e d[E] A dM, 

where Aj, i = 1, ... ,6 are some functions of the variables R, q, [E], and M. It is clear, 
therefore, that 

d9[ Tl ] A dQ [Tl ] A dQ [Tl] = -6Aid[E] A dT\ A dM A dT A dR A dg, 

and the points where the six-form vanishes coincide with those where A\ vanishes. A\ can 
be calculated from Eq. (fTJ) with the result 

Al= IVF + M^-4^- (45) 



Thus, A\ vanishes for 



IF] 2 - M4 a 2 



Do these points lie in the physical region that is given by inequality Q38|)? Let us substitute 
the value of [E] 2 into Eq. fl38|) and rewrite the result in the form 

M 2 , A g 2 \ 3M 2 
< 4 1 + 1 



i2 2 ~ I M 2 / M 2 + g 2 



This inequality implies inequality fl55|); hence, all zero points of Ai lie in the physical region. 
We conclude that [Ti] = const is no transversal surface. 

Let us try to shift [Ti] by a function dependent on the variables R, q, [E] and M: 

[T 1 ]=T 2 + X(R,q, [E],M). 

Our motivation for choosing such a form is that T 2 is then transformed by the A-symmetry 
fl55|) in the same way as [Ti] is, and that it is also a cyclic coordinate. Substitution for [Ti] 
into dO changes the coefficient A\ at di? A dg in d0T 2 by 

1 1 dq dR dR dq 
Comparing A±, Eq. d45|) , with dE /dR calculated from Eq. (p9j ), 

<9£ 1 , g 2 [£] 2 M 2 



dR 2 2M 2 2M 2 8Ti 2 ' 



we observe that 



dE q 2 

1 + m + Si*' 



which is always positive; thus, one possible shift is 
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[li] = T 2 - 2q. (46) 
If we perform this transformation in 0, we find out easily that another shift, 

H = T, + , (47) 

miraculously cancels some pesky cross terms in dO. Indeed, shifting the 'times' in the Cartan 
form, we have: 

-Ed[T x ) - [E\df x + MdJ = -EdT 2 - [E]dT 3 + MdJ 

_ ( B + g _ *hi _ ";u + *a r + _ « M . 

\ M 2 M 2 4RJ M 2 M 2 M 3 

Using Eqs. fl4"T| ) and (pB|), we obtain 

+ + m) dR+2 4( q+ St) dM - 

The last three terms do not contain other variables than q, R and M and so the terms with 
d[E] A dR, d[E] A dq and d[E] A dM disappear from dQ. We have not found any simple 
physical or geometrical interpretation of the new times T 2 and T3; observe that the shifts 
(||) and (g7D "mix the sides" . 

The fact that the last three terms in contain only three variables means that they can 
be simplified to just one term. For example, 



" q ;dq - ( q + ) dR + ^ ( q + 



2M 2 * V 2 M V M V 2M 

This suggests the second step of our transformation: 

p := ?M 2 + ^, (48) 

x := (49) 
We arrive so at the final shape of the Cartan form 

= -EdT 2 - [E]dT 3 + MdJ - pdx, (50) 



where E is given by Eq. ( PUD in which g and R are expressed in terms of p and x; such 
expressions can be obtained by solving Eqs. (|^) and (|49"D for g and i?. 

Consider Eq. (^). If we keep M fixed, p is an increasing function of q in the whole 
interval (—0x3,0x3), for 

^ = M 2 + 2q 2 > 0. 



Hence, the function maps the range (—00, 00) of q onto the range (— 00, 00) of p in a bijective 
way, and it possesses a unique inverse. One can explicitly write down this inverse using 
second and third roots; this will make E a complicated function of p, x, M and [E] (we shall 
write down this function later, cf. Eq. (|63D ). 
Eq. (pOD implies that 

dQ T2 A dQ Ti A dQ T2 = 6d[E] A dT 3 A dM A dT A dp A dar, 

rfe T3 A rfe T3 A rfe T3 = 6—— -d[E] A rfT 2 AdMAdJ AdpA dx, 

o[E\_ 

dE 

dQ T A dQ T A dQ T = -67^ d[E] A dT 3 A dM A dT 2 AdpA dx, 



and Eq. (p9|) gives 



d[E) M 2 
d£ M 

Thus, the surface T 3 = const is not transversal at [E] = 0, but T 2 = const and T = const 
are transversal everywhere. 

The Cartan form fl50|) is not yet very simple, because E is a complicated function of p, 
x, [E] and M. What we can still do is to extend the phase space to an eight-dimensional 
manifold with the cordinates E, [T 2 ], [E], T 2 , M, T, p and x, and to express the constraint 
between these variables as a polynomial. Indeed, applying Eq. fl49D , we can rewrite Eq. (|29|) 
as follows 

<f = A + — + [£] 2 "M 2 ; 
4x z x 

squaring Eq. (f48|) and substituting the above expression for q 2 into the result, we obtain 

" 2 = 5 (i? + T + ™' - M ' 2 ) (2? + T + ^ + M ' 2 ) 2 ' (51) 

This is a constraint that zs polynomial in the momenta p, E, [E] and M. The inequalities 
defining the ranges of the variables E, [T 2 ], [E], T 2 , M, T, p and x in the physical region are 
also simplified: 

[E] 2 < AE 2 , M > 0, x > 0. (52) 

Let us, finally, find the transformation of the variables p, x, M, T, E, T 2 , [E] and T 3 by 
the symmetry group. The results are as follows. We have for the A-transformation, 

p h-> p, x h-> MhM, ThT, E\-+E, [E)\-^[E], 

T2h " T2 "M' T3 ^ T3 "9M' T ^ T+ aM ; 
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for the time reflection a±, 

p^-p, x^x, M ^ M, Tk-T, E^E, [E}^[E], T 2 ^-T 2 , T 3 ^ -T 3 ; 
and for the left-right reflection a 2 , 

p^p, x^x, M i-> M, T h-> T, E^ E, [E] i-> -[£], T 2 ^ T 2 , T 3 ^ -T 3 . 



Thus, the simple form of the symmetry transformations is preserved by the shifts ( pESp and 
([17]) as well as by transformations (|48| ) and (f49[) . Eq. (|50|) implies that 

_ _ , fdA =, dA rrnl 9A r A \ 
by a A-transformation, so our extended A-transformation is a symmetry, 

e i-f -e 

by the time reflection o"i^| and 

0^0 

by the left-right reflection cr 2 . The constraint equation (KTI) is clearly invariant of the group. 
Thus, the action of the symmetry group is well-defined in the whole phase space and it 
preserves the physical region. 

VII. THE CONSTRAINT SURFACE 

In this section, we are going to study the properties of the constraint (|5lD. We are 
interested in the topology of the constraint surface; we shall investigate the salient question 
of whether or not any spurious solutions have been unintentionally included during the 
process of transforming the constraint to the polynomial form; and we shall look for solutions 



of the constraint (pi ) with respect to some momenta. 



Our new variational principle reads 

S = J dt(-ET 2 - [E]f 3 + MT-px- vxCi), (53) 



where 



ei=9pM^ + ^+[£] 2 -M 2 )(^ + ^ + 2[^ + M>) ■ (54) 



2 The sign change of G by o\ is necessary, because time reflections are anti-symplectic transform- 
ations. 
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The constraint surface is trivial in the direction of the three times T2, T3 and T, so we can 
limit ourselves to its projection to the five-dimensional space spanned by E, [E], M, p and 
x. Let us observe that C\ depends on E, [E] and x through a function that we shall call B, 

and that the gradient of B is non-zero. It is, therefore, advantageous, to write C\ as a 
composed function, 

d = 9p 2 -{B- M 2 )(25 + M 2 ) 2 . 

Let us first look for the singular points of the constraint surface (that is, where the 
gradient of C\ vanishes). We have 

f = 18 "- < 56 > 

1 6M 3 (25+M 2 ), (57) 



<9M 
dCi 
~dB 



-3(25 + M 2 )(2B- M 2 ). (58) 



As all three derivatives must vanish at a singular point, all such points are determined by 
the equations 

p = 0, 25 + M 2 = 0. (59) 

All points satisfying Eqs. ( |55D are solutions of Eq. (CT); however, the physical ranges ( |5"2"D of 
the variables E 1 , [£7] and 2 do not allow B to be negative. We conclude that the constraint 
surface is regular (smooth) in the physical region fl5"2"|). 

The topology of the constraint surface can be found in the shortest way, if we consider 
C\ as a function Ci(B) of B keeping all other variables constant. Studying Eqs. (|54] ) and 
(|58|) we can easily see that C\{B) decreases in the interval (—00, — M 2 /2) from 00 to 9p 2 , it 
increases in the interval (— M 2 /2), M 2 /2), from 9p 2 to 9p 2 + 2M 6 and it decreases again in the 
interval (M 2 /2, 00) from 9p 2 + 2M 6 to —00. It follows that there is only one solution for p ^ 
(because C X {B) = 9p 2 for B = M 2 ). As Ci(M 2 ) = 9p 2 > 0, the solution must satisfy the 
inequality B > M 2 . This solution depends smoothly on the values of all other variables. For 
p = 0, there are, however, two solutions: one with B = M 2 , which is a continuous extension 
of the previous p 7^ case, and a new one with B = — M 2 /2, which appears when the local 
minimum at — M 2 /2 touches the I?-axis. This shows that the constraint surface consists of 
two components: one is a seven-dimensional smooth manifold, let us denote it by T 7 , that 
lies in the region B > M 2 , and one is a six-dimensional smooth manifold, let us denote it by 
r 6 , that is defined by Eqs. (|59|) . 

Next, let us study the solvability of the constraint equation ([H]) with respect to the 
conserved momenta. The unique solution of the constraint (]51f) with respect to B at T 7 is 
given by 

B = f(p,M), (60) 
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where 



/(p, M) := ^9p 2 + M 6 + 3pyV + 2M^ + ~\/9p 2 + M 6 - + 2M 6 . (61) 



The equation ( pT|) has, in fact, three independent solutions, but only one of them is real. 
For p = 0, B reaches its minimum M 2 on T 7 ; thus, at T-j, 

B > M 2 . (62) 

Given B, one can solve uniquely for E and there is also a unique solution for [E] 2 . The 
solution for E is 

E = _J__x\E? + ^/(p,M). (63) 

Similarly, we can find the solution of the constraint ([y]) with respect to M. Consider C\ 
as a function Ci(M) of M keeping all other variables fixed. Ci(M) is symmetric with respect 
to the reflection on the Ci-axis. We have to distinguish two cases: B < 0, B = 0: Ci(M) 
decreases in the interval (— oo, 2B) from oo to 9p 2 , it increases in the interval (2B, 0) from 
9p 2 to the local maximum 9p 2 — AB 3 at M = 0, it decreases again in the interval (0, —2B) 
from 9p 2 — AB 3 to 9p 2 , and, finally, it increases in (— 2B, oo) from 9p 2 to oo. A solution 
exists only for p = 0, and there are two solutions, then, M = ±V — 2B. This is at the surface 
r 6 . B > : then, there is only one minimum, at M = 0, with the value 9p 2 — AB 3 ; there are 
no other extrema. Consider Eq. ([Hp as an equation for M~ 2 . There are three real solutions 
for M -2 , but all of them are negative unless 

AB 3 - 9p 2 > 0, (64) 

and then there is only one non-negative solution. A reasonable solution is, therefore, unique, 
it depends continuously of the other variables and it lies at T 7 . We can write it in the form 

M = ± Jjg^/ 1 _ 3 P V 1 ' 2 



AB V 3 - 9p 2 ) 

This solution can, of course, be also expressed by means of square and third roots, but only 
if one employs complex numbers. We can see that squaring of M introduced non-physical 
solutions with negative M; we have to choose the positive branch. Observe that all points 
at the surface T 7 must satisfy the inequalities fl5"^ ) and (Q). 

The equations of motion that follow from the action (^) can be written as follows 

(65) 
(66) 

(67) 

(68) 
(69) 
(70) 



E = 


[E] — M — 0. 


f 2 = 


dCi dB 


Vx dB BE' 


f 3 = 


dC x dB 


Ul dB d[E] 


t = 


dd 






dCx dB 


p = 


1 dB dx ' 




dd 


x = 


dp 



22 



It is easy to derive the radial equation from them. On Tq, we have only a static solutions: 



x = T = ft = T 2 = T 3 = 0. 
On r 7 , t > for v x > and Eqs. Q57|), and yield 

dT = _ M 3 (2B + M 2 ) ' 



Using the constraint (^l|) to exclude p, we obtain the radial equation on IV 



Tf-^wiw- 1 - (71) 



Let us compare it with Eq. dl3|) . Eqs. d65|) and (|49|) imply that 



whereas Eqs. fl55|), (||) and ( P^ ) give 



1 c?i? 



Hence, the radial equations (|71|) and (|X3|) are equivalent, and the dynamics of the new 
action on coincides with that of the old action. T 6 consists of the unintentionally added 
unphysical solutions, at least on the phase space defined by Eqs. dl8|), but nothing is added 
in the physical region. 

Finally, we shall study the question of the monotonicity of the time functions T 2 , T3 and 



T along the dynamical trajectories at T 7 . Eqs. (|57|) and (p8|) show that the right-hand sides 



of Eqs. (^Bp and flog ) cannot vanish, and we have for V\ > 0: 

T 2 > 0, t > 0. 

Hence, T 2 and T are good times and, as we have seen above, the constraint can be solved 
uniquelly for the corresponding momenta, E and M. As for T3, Eq. ( |67D implies that T3 > 
for [E] > 0, T 3 = for [E\ = and T 3 < for [E\ < 0. The solvability with respect to [E\ 
is also only partial: there is a unique positive and a unique negative solution. Thus T 3 is a 
good time only in some special the next section will show. 



VIII. SOME INTERESTING SPECIAL CASES 

In this section, we shall remove some degrees of freedom and describe the resulting 
simpler models in terms of the new variables. 
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A. Dust degrees of freedom removed 



Here, we remove the variables M and T and return to the system considered in Ref. [TH 



It has two degrees of freedom: the black hole mass E_ and the position x of the shell. To 
this aim, we first choose a particular value of M and demote it so to a mere parameter. 
In this way, a submanifold T 7 m of T 7 emerges. Second, we take a quotient of T 7 m by the 
T-curves (see Ref. [[y| for more detail of this exclusion of a cyclic variable). The system 
that results has the action 



Sm = J dt(-ET 2 - [E]f 3 - px - v 1 C 1 



where C\ is given by Eq. (p4|) as before, only M is a parameter now. The equations of motion 
comprise Eq. (@) without M = 0, as well as Eqs. Q6?j), fl67|), © and (^). 

The properties of this system are analogous to that of the original one. In particular, the 
constraint is regular everywhere on T 7 m because the derivative (^) and ( |55| ) cannot both 
simultaneously vanish, T 2 is a good time variable and the constraint equation is uniquely 
solvable for E. 



B. Black hole degree of freedom removed 

Here, we demote the variable E_ to a parameter and remove the corresponding cyclic 
variable T 2 /2 — T 3 ; the remaining two degrees of freedom are the internal energy M and the 
position x of the shell. The system obtained in this way will describe the dynamics of the 
dust shell in the field of a fixed black hole if EL > or in the Minkowski space if E_ = 0. 
One returns so to the system considered in Ref. 
has been performed in Ref. |L5|. 



13]; the procedure, in different coordinates, 



To begin with, we transform variables as follows: 



E 



E+ + E. 



1 2 



T 2 /2 + T 3 , 



[E] = E + — E_, T 2 -=T 2 /2-T 3 . 

Then, we choose a particular value for E_; this defines the black hole mass, and formally, a 
six-dimensional submanifold T 7 e of T 7 . Finally, we take the quotient T 7 e/T_ of Y 7 e by the 
T_-curves. The result is the action 

S E = J dt(-E + f 2 + + MT - px - v x C 2 ), 

where 

C 2 = 9p 2 -( J B 2 -M 2 )(2 J B 2 + M 2 ), 

and 
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B -i+^+< £ +- £ -> 2 

The variable T is a good time and C% = is of course solvable for M exactly as in the 
general case. The variable T 2 + is, for large values of the parameter E_, not a good time, 
because it is a combination T 2 /2 + T 3 of T 2 that increases along all dynamical trajectories 
and T 3 that increases along those with E + > E_ and decreases along those with E + < E_. 
One can show that T 2 + is positive for E + — E_ > and negative for E + — E_ < — M and that 
there is no surface for E_ > M that would be transversal to both the dynamical trajectories 
and the orbits of the continuous subgroup. f\ The derivative of C 2 with respect to E + changes 
sign somewhere if E_ is sufficiently large and so there are more that one, or no solution of 
the constraint with respect to E + . 

There is one prominent exception: suppose that E- = 0. Then, E + — E_ is positive if 
E + is. Hence, in the physical region (where E + > 0), T 2 + is a good time and the constraint 
has a unique solution for E + there, namely 

E + = ~ + sfttp, M), 

where / is defined by Eq. (|6lD. 

One can also remove all four variables T, M, E- and T 2 _ so that only one degree of 
freedom, the position x of the shell remains, but it is not difficult to work out the properties 
of the system using the results obtained in this section, and we left it as an easy exercise to 
the reader. 

The special case E_ = is also interesting, because it is a part of the boundary of the 
physical region given by inequalities (p2[) on Fj. The variational principle (|53|) does not 
break down at the boundary, that is, it defines a regular dynamics there. From this, it must 
clearly follow that the action defines a regular dynamics even at those points of I\ that do 
not satisfy the conditions (|52|). We have studied the meaning of this dynamics and of these 
points in Sec. [TJ. 
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3 Still, one can transform T 2 + so that the resulting time levels are transversal to the dynamical 
trajectories, but a mere shift is not sufficient; one must screw the level to a helix form with the 
helix axis in the -E^-direction; this may be done by a singular "shift" proportional to {E + — E^)^ 1 . 
We shall not go into detail. 
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